What Is Claimed Is: 



1 . A method for using a computer system to solve a global inequality 
constrained optimization problem specified by a function / and a set of inequality 
constraints p^x) < 0 (i=l, wherein / and p t are scalar functions of a vector 
x = (xj, X2, ... x„) 9 the method comprising: 

receiving a representation of the function / and the set of inequality 
constraints at the computer system; 

storing the representation in a memory within the computer system; 

performing an interval inequality constrained global optimization process 
to compute guaranteed bounds on a globally minimum value of the function/x) 
subject to the set of inequality constraints; 

wherein performing the interval inequality constrained global optimization 
process involves, 

applying term consistency to a set of relations associated 
with the global inequality constrained optimization problem over a 
subbox X, and excluding any portion of the subbox X that violates 
any of these relations, 

applying box consistency to the set of relations associated 
with the global inequality constrained optimization problem over 
the subbox X, and excluding any portion of the subbox X that 
violates any of these relations, and 

performing an interval Newton step for the global 
inequality constrained optimization problem over the subbox X to 
produce a resulting subbox Y, wherein the point of expansion of 
the interval Newton step is a point x. 
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2. The method of claim 1, wherein applying term consistency to the 
set of relations involves applying term consistency to the set of inequality 
constraints p { (x) < 0 (/=/, ...,m) over the subbox X. 

3. The method of claim 1, wherein applying box consistency to the 
set of relations involves applying box consistency to the set of inequality 
constraints p t {x) < 0 (i=7, ...,m) over the subbox X. 

4. The method of claim 1 , 

wherein performing the interval inequality constrained global optimization 
process involves, 

keeping track of a smallest upper bound fjbar of the 
function f(x) at a feasible point x, 

removing from consideration any subbox X for which 
/X) >f_ban 

wherein applying term consistency to the set of relations involves applying 
term consistency to the f_bar inequality fix) < fjbar over the subbox X. 

5. The method of claim 4, wherein applying box consistency to the 
set of relations involves applying box consistency to the fjbar inequality 

fix) < fjbar over the subbox X. 

6. The method of claim 1, wherein if the subbox X is strictly feasible 
(p t (X) < 0 for all ...,n), performing the interval inequality constrained global 
optimization process involves: 

determining a gradient g(x) of the function/(x), wherein g(x) includes 
components g t (x) (i=l, ...,/?); 
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removing from consideration any subbox for which g(x) is bounded away 
from zero, thereby indicating that the subbox does not include an extremum of 
J{x); and 

wherein applying term consistency to the set of relations involves applying 
term consistency to each component g,(x)=0 (z=l, ...,ri) of g(x)=0 over the subbox 
X. 

7. The method of claim 6, wherein applying box consistency to the 
set of relations involves applying box consistency to each component 
gi(x)=0 (/=/, ...,«) of g(x)=0 over the subbox X. 

8. The method of claim 1 , wherein if the subbox X is strictly feasible 
(Pi(X) < 0 for all i=l, ...,/?), performing the interval inequality constrained global 
optimization process involves: 

determining diagonal elements H u (x) (i=l,...,ri) of the Hessian of the 
function^*); 

removing from consideration any subbox for which a diagonal element 
H n (X) of the Hessian over the subbox X is always negative, indicating that the 
function /is not convex over the subbox X and consequently does not contain a 
global minimum within the subbox X; and 

wherein applying term consistency to the set of relations involves applying 
term consistency to each inequality H n (x) > 0 (7=1, ...,ri) over the subbox X. 

9. The method of claim 8, wherein applying box consistency to the 
set of relations involves applying box consistency to each inequality 

H H (x) ^ 0 (7=1, ...,ri) over the subbox X. 
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1 0. The method of claim 1 , wherein if the subbox X is strictly feasible 
(p,(X) < 0 for all f= 7, .„,«), performing the interval Newton step involves: 

computing the Jacobian J(x,X) of the gradient of the function/evaluated 
with respect to a point x over the subbox X; and 

computing an approximate inverse B of the center of J(x,X), 
using the approximate inverse B to analytically determine the system Bg(x), 
wherein g(x) is the gradient of the function/x), and wherein g(x) includes 
components g t (x) (7=1, ...,«). 

1 1 . The method of claim 1 0, wherein applying term consistency to the 
set of relations involves applying term consistency to each component 
(Bg(x)), = 0 (z"=7, ...,«) to solve for the variable x l over the subbox X. 

12. The method of claim 10, wherein applying box consistency to the 
set of relations involves applying box consistency to each component 
(Bg(x))/ = 0 (7=7, ...,ri) to solve for the variable x ; - over the subbox X. 

13. The method of claim 1 , wherein performing the interval Newton 
step involves performing the Newton step on the John conditions. 

14. The method of claim 1 , 

wherein performing the interval inequality constrained global optimization 
process involves, 

linearizing the set of inequality constraints to produce a set 
of linear inequality constraints with interval coefficients that 
enclose the nonlinear inequality constraints, and 
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1 preconditioning the set of linear inequality constraints 

2 through additive linear combinations to produce a set of 

3 preconditioned linear inequality constraints; and 

4 wherein applying term consistency to the set of relations involves applying 

5 term consistency to the set of preconditioned linear inequality constraints over the 

6 subbox X. 



1 15. The method of claim 1 4, wherein applying box consistency to the 

2 set of relations involves applying box consistency to the set of preconditioned 

3 linear inequality constraints over the subbox X. 



1 16. The method of claim 1 , wherein applying term consistency 

2 involves: 

3 symbolically manipulating an equation within the computer system to 

4 solve for a term, g(x }), thereby producing a modified equation g(x )) = h{\\ 

5 wherein the term g(x }) can be analytically inverted to produce an inverse function 

6 g\y); 

7 substituting the subbox X into the modified equation to produce the 

8 equation g(X = h(X) ; 

9 solving for X ) = g-\h(X)); and 

1 0 intersecting X ) with the >th element of the subbox X to produce a new 

1 1 subbox X + ; 

12 wherein the new subbox X + contains all solutions of the equation within 

13 the subbox X, and wherein the size of the new subbox X + is less than or equal to 

14 the size of the subbox X. 
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1 17. The method of claim 1 , wherein performing the interval Newton 

2 step involves: 

3 computing J(x,X), wherein J(x,X) is the Jacobian of the function f 

4 evaluated as a function of x over the subbox X; and 

5 determining if J(x,X) is regular as a byproduct of solving for the subbox Y 

6 that contains values of y that satisfy M(x,X)(y-x) = r(x), where 

7 M(x ? X) = BJ(x,X), r(x) = -Bf(x), and B is an approximate inverse of the center of 

8 J(x,X). 



1 1 8. A computer-readable storage medium storing instructions that 

2 when executed by a computer cause the computer to perform a method for using a 

3 computer system to solve a global inequality constrained optimization problem 

4 specified by a function / and a set of inequality constraints pfa) < 0 (1=1, 

5 wherein / and p t are scalar functions of a vector x = (xj, x% x 3> ... x n \ the method 

6 comprising: 

7 receiving a representation of the function / and the set of inequality 

8 constraints at the computer system; 

9 storing the representation in a memory within the computer system; 

10 performing an interval inequality constrained global optimization process 

1 1 to compute guaranteed bounds on a globally minimum value of the function/(x) 

12 subject to the set of inequality constraints; 

13 wherein performing the interval inequality constrained global optimization 

14 process involves, 

15 applying term consistency to a set of relations associated 

16 with the global inequality constrained optimization problem over a 

17 subbox X 5 and excluding any portion of the subbox X that violates 

1 8 any of these relations, 
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1 9 applying box consistency to the set of relations associated 

20 with the global inequality constrained optimization problem over 

21 the subbox X, and excluding any portion of the subbox X that 

22 violates any of these relations, and 

23 performing an interval Newton step for the global 

24 inequality constrained optimization problem over the subbox X to 

25 produce a resulting subbox Y, wherein the point of expansion of 

26 the interval Newton step is a point x. 



1 19. The computer-readable storage medium of claim 1 8, wherein 

2 applying term consistency to the set of relations involves applying term 

3 consistency to the set of inequality constraints p l (x) < 0 (i=l, ...,m) over the 

4 subbox X. 



1 20. The computer-readable storage medium of claim 1 8, wherein 

2 applying box consistency to the set of relations involves applying box consistency 

3 to the set of inequality constraints p{x) < 0 (zW, ...,m) over the subbox X. 



1 21. The computer-readable storage medium of claim 1 8, 

2 wherein performing the interval inequality constrained global optimization 

3 process involves, 

4 keeping track of a smallest upper bound fjbar of the 

5 function/x) at a feasible point x, 

6 removing from consideration any subbox X for which 

7 /X) >fjar; 

8 wherein applying term consistency to the set of relations involves applying 

9 term consistency to the fjbar inequality fix) < f_bar over the subbox X. 
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1 22. The computer-readable storage medium of claim 21 , wherein 

2 applying box consistency to the set of relations involves applying box consistency 

3 to the fjbar inequality fix) < f_bar over the subbox X. 



1 23 . The computer-readable storage medium of claim 22 , wherein if the 

2 subbox X is strictly feasible 0/(X) < 0 for all i=l, ...,n), performing the interval 

3 inequality constrained global optimization process involves: 

4 determining a gradient g(x) of the function/x), wherein g(x) includes 

IL.~. 

hi 5 components g,(x) (*'=i, ...,«); 

^ 6 removing from consideration any subbox for which g(x) is bounded away 

_ » ^ ^ 

fy 7 from zero, thereby indicating that the subbox does not include an extremum of 

8 J(x); and 

9 wherein applying term consistency to the set of relations involves applying 

1 0 term consistency to each component gt(x)=0 (i= 1, . . n) of g(x)=0 over the subbox 

11 X. 



1 24. The computer-readable storage medium of claim 23 , wherein 

2 applying box consistency to the set of relations involves applying box consistency 

3 to each component g,(x>K) (i=l, ...,«) of g(x)=0 over the subbox X. 



1 25 . The computer-readable storage medium of claim 1 8 ? wherein if the 

2 subbox X is strictly feasible (p^X) < 0 for all ...,ri), performing the interval 

3 inequality constrained global optimization process involves: 

4 determining diagonal elements H h (x) (7=1, of the Hessian of the 

5 function fix); 



35 

Attorney Docket No. SUN-P6432-SPL Inventors: Walster et al. 

ARPHASUN MICROSYSTEMS\SUN-P6432-SPL\SUN-P6432-SPL APPLICATION DOC 



6 removing from consideration any subbox for which a diagonal element 

7 Hu(X) of the Hessian over the subbox X is always negative, indicating that the 

8 function /is not convex over the subbox X and consequently does not contain a 

9 global minimum within the subbox X; and 

1 0 wherein applying term consistency to the set of relations involves applying 

1 1 term consistency to each inequality H u (x) > 0 (7=1, ...,ri) over the subbox X. 

1 26. The computer-readable storage medium of claim 25 , wherein 

2 applying box consistency to the set of relations involves applying box consistency 

3 to each inequality #„(x) > 0 (7=1, ...,«) over the subbox X. 

1 27. The computer-readable storage medium of claim 1 8, wherein if the 



2 subbox X is strictly feasible (p,(X) < 0 for all i=l, performing the interval 

3 Newton step involves: 

4 computing the Jacobian J(x,X) of the gradient of the function / evaluated 

5 with respect to a point x over the subbox X; and 

6 computing an approximate inverse B of the center of J(x,X), 

7 using the approximate inverse B to analytically determine the system Bg(x), 

8 wherein g(x) is the gradient of the function/x), and wherein g(x) includes 

9 components g,(x) (i = 1, . . . , n). 

1 28. The computer-readable storage medium of claim 27, wherein 

2 applying term consistency to the set of relations involves applying term 

3 consistency to each component (Bg(x)), = 0 (7=7, ...,n) to solve for the variable x t 

4 over the subbox X. 
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1 29. The computer-readable storage medium of claim 27, wherein 

2 applying box consistency to the set of relations involves applying box consistency 

3 to each component (Bg(x)), = 0 (7= 7, ...,«) to solve for the variable x t over the 

4 subbox X. 



1 30. The computer-readable storage medium of claim 1 8, wherein 

2 performing the interval Newton step involves performing the interval Newton step 

-> 

3 on the John conditions. 



1 31. The computer-readable storage medium of claim 1 8, 

2 wherein performing the interval inequality constrained global optimization 

3 process involves, 

4 linearizing the set of inequality constraints to produce a set 

5 of linear inequality constraints with interval coefficients that 

6 enclose the nonlinear inequality constraints, and 

7 preconditioning the set of linear inequality constraints 

8 through additive linear combinations to produce a set of 

9 preconditioned linear inequality constraints; and 

1 0 wherein applying term consistency to the set of relations involves applying 

1 1 term consistency to the set of preconditioned linear inequality constraints over the 

1 2 subbox X. 



1 32. The computer-readable storage medium of claim 3 1 , wherein 

2 applying box consistency to the set of relations involves applying box consistency 

3 to the set of preconditioned linear inequality constraints over the subbox X. 
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1 33. The computer-readable storage medium of claim 1 8 ? wherein 

2 applying term consistency involves: 

3 symbolically manipulating an equation within the computer system to 

4 solve for a term, g(x )\ thereby producing a modified equation g(x )) = h(x), 

5 wherein the term g(x }) can be analytically inverted to produce an inverse function 

6 g'(y); 

7 substituting the subbox X into the modified equation to produce the 

8 equation g(X )) = h(X); 

9 solving for X ) = g _/ (/z(X)); and 

1 0 intersecting X ) with the y-th element of the subbox X to produce a new 

1 1 subbox X + ; 

£ 12 wherein the new subbox X + contains all solutions of the equation within 

JJj 1 3 the subbox X, and wherein the size of the new subbox X is less than or equal to 

m 14 the size of the subbox X. 



1 34. The computer-readable storage medium of claim 1 8, wherein 

2 performing the interval Newton step involves: 

3 computing J(x,X) 5 wherein J(x,X) is the Jacobian of the function f 

4 evaluated as a function of x over the subbox X; and 

5 determining if J(x,X) is regular as a byproduct of solving for the subbox Y 

6 that contains values of y that satisfy M(x,X)(y-x) = r(x), where 

7 M(x 5 X) = B J(x 5 X), r(x) = -Bf(x), and B is an approximate inverse of the center of 

8 J(x,X). 



1 35 . An apparatus that solves a global inequality constrained 

2 optimization problem specified by a function / and a set of inequality constraints 
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3 jt?i(x) < 0 (i=l, ..,4 wherein / and p z are scalar functions of a vector 

4 x = (x/, x 2y x 3f ... x„), the apparatus comprising: 

5 a receiving mechanism that is configured to receive a representation of the 

6 function/ and the set of inequality constraints at the computer system; 

7 a memory for storing the representation; 

8 an interval global optimization mechanism that is configured to perform 

9 an interval inequality constrained global optimization process to compute 

10 guaranteed bounds on a globally minimum value of the function/(x) subject to the 

1 1 set of inequality constraints; 

12 a term consistency mechanism within the interval global optimization 

13 mechanism that is configured to apply term consistency to a set of relations 

14 associated with the global inequality constrained optimization problem over a 

1 5 subbox X ? and to exclude any portion of the subbox X that violates any of these 

1 6 relations, 

17 a box consistency mechanism within the interval global optimization 

1 8 mechanism that is configured to apply box consistency to the set of relations 

19 associated with the global inequality constrained optimization problem over the 

20 subbox X, and to exclude any portion of the subbox X that violates any of these 

21 relations, and 

22 an interval Newton mechanism within the interval global optimization 

23 mechanism that is configured to perform an interval Newton step for the global 

24 inequality constrained optimization problem over the subbox X to produce a 

25 resulting subbox Y, wherein the point of expansion of the interval Newton step is 

26 a point x. 
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36. The apparatus of claim 35, wherein the term consistency 
mechanism is configured to apply term consistency to the set of inequality 
constraints p^x) < 0 (i=l, ...,m) over the subbox X. 

37. The apparatus of claim 35, wherein the box consistency 
mechanism is configured to apply box consistency to the set of inequality 
constraints p,( x ) * 0 (i=l, ...,m) over the subbox X. 

38. The apparatus of claim 35, 

wherein the interval global optimization mechanism is configured to, 

keep track of a smallest upper bound f_bar of the function 
fix) at a feasible point x, and to 

remove from consideration any subbox X for which 

/X) >fjar; 

wherein the term consistency mechanism is configured to apply term 
consistency to the f_bar inequality fix) <f_bar over the subbox X. 

39. The apparatus of claim 38, wherein the box consistency 
mechanism is configured to apply box consistency to the f_bar inequality 
fix) < fjbar over the subbox X. 

40. The apparatus of claim 35, wherein if the subbox X is strictly 
feasible (p,{X) < 0 for all the interval global optimization mechanism is 
configured to: 

determine a gradient g(x) of the function/x), wherein g(x) includes 
components g^x) ( z== ^ — > n )l 
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1 remove from consideration any subbox for which g(x) is bounded away 

2 from zero, thereby indicating that the subbox does not include an extremum of 

3 f(x); and 

4 the term consistency mechanism is configured to apply term consistency to 

5 each component gi(x)=0 (i=l, ...,ri) of g(x)=0 over the subbox X. 



1 41 . The apparatus of claim 40, wherein the box consistency 

2 mechanism is configured to apply box consistency to each component 

3 giOO^O ...,n)of g(x)=0 over the subbox X. 



1 42. The apparatus of claim 35, wherein if the subbox X is strictly 

2 feasible (p,(X) < 0 for all i=l 9 ...,«), the interval global optimization mechanism is 

3 configured to: 

4 determine diagonal elements H n (x) (i=l, ...,ri) of the Hessian of the 

5 functionXx); 

6 remove from consideration any subbox for which a diagonal element 

7 H U (X) of the Hessian over the subbox X is always negative, indicating that the 

8 function /is not convex over the subbox X and consequently does not contain a 

9 global minimum within the subbox X; and 

10 the term consistency mechanism is configured to apply term consistency to 

1 1 each inequality H u (x) z 0 (z=l, over the subbox X. 



1 43 . The apparatus of claim 42, wherein the box consistency 

2 mechanism is configured to apply box consistency to each inequality 

3 H u (x) > 0 (/==1, ...,ri) over the subbox X. 



41 



Attorney Docket No. SUN-P6432-SPL Inventors: Walster et al. 

ARPH:\SUN MICR0SYSTEMS\SUN-P6432-SPL\SUN-P6432-SPL APPLICATION2.DOC 



44. The apparatus of claim 35, wherein if the subbox X is strictly 
feasible (p,(X) < 0 for all /=/, the interval global optimization mechanism is 
configured to perform the interval Newton step by: 

computing the Jacobian J(x,X) of the gradient of the function / evaluated 
with respect to a point x over the subbox X; and 

computing an approximate inverse B of the center of J(x ? X), 
using the approximate inverse B to analytically determine the system Bg(x), 
wherein g(x) is the gradient of the function/x), and wherein g(x) includes 
components g t (x) (z=l, 

45. The apparatus of claim 44, the term consistency mechanism is 
configured to apply term consistency to each component (Bg(x)) 7 = 0 (1=1, ...,ri) to 
solve for the variable x, over the subbox X. 

46. The apparatus of claim 44, the box consistency mechanism is 
configured to apply box consistency to each component (Bg(x)),- = 0 (i=l, ...,«) to 
solve for the variable x, over the subbox X. 

47. The apparatus of claim 35, wherein the interval Newton 
mechanism is configured to perform the Newton step on the John conditions. 

48. The apparatus of claim 35, 

wherein the interval global optimization mechanism is configured to: 

linearize the set of inequality constraints to produce a set of 
linear inequality constraints with interval coefficients that enclose 
the nonlinear inequality constraints, and to 
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1 precondition the set of linear inequality constraints through 

2 additive linear combinations to produce a set of preconditioned 

3 linear inequality constraints; and 

4 wherein the term consistency mechanism is configured to apply term 

5 consistency to the set of preconditioned linear inequality constraints over the 

6 subbox X. 



1 49. The apparatus of claim 48, wherein the box consistency 

2 mechanism is configured to apply box consistency to the set of preconditioned 

3 linear inequality constraints over the subbox X. 



1 50. The apparatus of claim 35, wherein the term consistency 

2 mechanism is configured to : 

3 symbolically manipulate an equation within the computer system to solve 

4 for a term, g(x thereby producing a modified equation g(x }) = h(x), wherein 

5 the term g(x }) can be analytically inverted to produce an inverse function g'\y); 

6 substitute the subbox X into the modified equation to produce the equation 

7 g(X>j) = h(X); 

8 solve for X ) = g ~\h(X)); and 

9 intersect X ) with the y'-th element of the subbox X to produce a new 

10 subbox X ; 

1 1 wherein the new subbox X + contains all solutions of the equation within 

12 the subbox X, and wherein the size of the new subbox X + is less than or equal to 

1 3 the size of the subbox X. 



1 51. The apparatus of claim 35, wherein the interval Newton 

2 mechanism is configured to: 
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compute J(x,X), wherein J(x,X) is the Jacobian of the function f evaluated 
as a function of x over the subbox X; and to 

determine if J(x,X) is regular as a byproduct of solving for the subbox Y 
that contains values of y that satisfy M(x,X)(y-x) = r(x), where 

M(x,X) = BJ(x,X), r(x) = -Bf(x), and B is an approximate inverse of the center of 
J(x,X). 
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